Engineering two-mode interactions in ion traps 
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We describe how two vibrational degrees of freedom of a single trapped ion can be coupled 
through the action of suitably-chosen laser excitation. We concentrate on a two-dimensional ion 
trap with dissimilar vibrational frequencies in the x- and j/-directions of motion, and derive from 
first principles a variety of quantized two-mode couplings, concentrating on a linear coupling which 
takes excitations from one mode to another. We demonstrate how this can result in a state rotation, 
in which it is possible to transfer the motional state of the ion from say the a;-direction to the 
y-direction without prior knowledge of that motional state. 
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I. INTRODUCTION 

In recent years advances in the cooling and trapping of 
ions have led to a situation in which the centre-of-mass 
(CM) motion of trapped ions has to be treated quan- 
tum mechanically |^]. This motion can be coherently 
controlled by coupling the ion's external and internal de- 
grees of freedom through laser irradiation |^-^. Sys- 
tems of trapped ions have been employed to demonstrate 
experimentally the generation and measurement of non- 
classical states of the ion 's CM motion [^^. Further- 
more trapped ions have been used to implement quantum 



logic gates 11-13 



Most of the previous investigations have focussed on 
the one-dimensional quantum motion of trapped ions. 
Recently, Gou et. al. JT^-p^ considered the generation 
of particular two-mode states of an ion. In this paper we 
address the issue of how to engineer a class of interac- 
tions between two of the quantized motional degrees of 
freedom of a single trapped ion. We assume that the ion 
is confined within a trap potential that can be closely 
approximated by a two-dimensional harmonic well. In 
this case the CM motion of the ion is completely equiv- 
alent to that of a two-dimensional harmonic oscillator, 
characterized by two frequencies of oscillation i>a and i>h 
in orthogonal directions x and y, and the corresponding 
operators al (a) and w (b) create (annihilate) vibrational 
excitations in the x- and y-directions. The interaction 
that we want to engineer is of the parametric form 
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(1) 



where fc^ and kb are positive integers and g is a complex 
coupling constant. In particular, we note that the pow- 
ers ka, and fcf,, can be independently controlled to take 
on any positive integer numbers, and the phase of the 
coupling constant is freely adjustable. To give specific 
examples of this class of interactions between the two vi- 
brational modes a and b, we address the two coupling 
Hamiltonians, 



(3) 



(2) 
(3) 



The Hamiltonian (g) generates an active rotation of the 
two-dimensional quantized motional state of the ion at 
a frequency g, where g is real. Here ifj^^ is the kind 
of Hamiltonian associated with a linear coupler or beam 
splitter in optics (see e.g. jl^ and refs. therein). There, 
a photon in mode a is annihilated and a photon in mode 
b is created, and vice versa. In a trapped ion, vibra- 
tional anticorrelated SU(2) states of motion characteris- 
tic of this kind of linear coupling can be generated . 

The linear coupling, makes it possible to transfer 

the motional state of the ion from, say the x-direction 
into the y-direction without prior knowledge of that mo- 
tional state and irrespective of whether it is a pure or a 
mixed state. In the situation in which one may want to 
use the quantized motion in the x-direction for quantum 
computation [ pd] , perhaps later entangling the quan- 
tum state of motion with internal electronic states, the 
y-direction can then be employed as a quantum mem- 
ory element. Note the key point here is that states of 
motion in the a;-direction can be transferred entirely to 

the y-direction without reading out their nature, entirely 

" (3) 

non-destructively. The Hamiltonian iJj is of the three- 
photon down-conversion kind: in optics, it represents a 
process in which one pump photon in mode b is annihi- 
lated and three photons in mode a are created, and vice 
versa. This process is known to be highly peculiar: un- 
like its two-photon down-conversion cousin, quantization 
of the pump is essential to avoid pathological divergences 
pO| . These are avoided in a fully quantized treatment, 
where the pump and down-converted field modes become 
highly entangled pl| . 

In section we first introduce a two-mode Raman 
transition which couples the electronic and motional de- 
grees of freedom of the ion. Choosing the initial state of 
the ion to be a direct product of an arbitrary motional 
state and a specific electronic state, we then decouple 
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the electronic and motional dynamics of the ion through 
a particular configuration of laser beams (section . In 
the Lamb-Dicke approximation and in the limit of suit- 
able trap anisotropy we obtain the above Hamiltonian (Q) 
for various sideband detunings of the lasers. We then ex- 
amine the severity of the approximations made to obtain 
the Hamiltonian (|l|). In section |^ we obtain analytical 
estimates regarding the effects of off-resonant and higher 
on-resonant processes. In section ^ we specialize to the 
case ka = kh = 1, and show that the Hamiltonian H^^"^ 
rotates the motional quantum state of the ion. Finally, 
we perform a numerical analysis of the complete quan- 
tum dynamics and find that the Hamiltonian (|^) can be 
accurately engineered over a range of parameters. 



II. GENERAL TWO-MODE RAMAN COUPLING 

In the following we describe the Raman coupling 
which we use to engineer the Hamiltonian given in 
(|I|). We consider an effective three- level ion in a A- 
configuration, confined within a two-dimensional har- 
monic trap as illustrated in Fig.|^. The trap is char- 
acterized by the two frequencies Va and Vb which de- 
scribe the harmonic potential in x- and y-direction re- 
spectively. As shown in Fig.|l|, the ion is excited by 
two linearly polarized laser beams which propagate in 
the X- and y-directions connecting levels |1) <^ |2) and 
1 2) |3). These beams are far detuned from the ex- 
cited state 1 2), in order to generate a stimulated Raman 
transition between the two states |1) and |3). We assume 
states |1) and |3) to be ground state hyperfine sublevels. 



We do not include decoherence effects in our model for 
the following reason. The Raman coupled energy level 
scheme greatly suppresses the spontaneous emission be- 
tween the two ground state levels |3) and |1) as these 
states are coupled by Ml and E2 transitions at best. At 
the same time we neglect the effects of spontaneous emis- 
sion from level |2), as the coupling to the excited state 
can be effectively eliminated over the timescales of in- 
terest to us here when the laser beams are far detuned. 
Another source of decoherence in ion trap experiments is 
classical noise in the laser beams and trapping potential. 
This may be described using so-called intrinsic decoher- 
ence models, e.g. |2^], of dephasing. The effects of this 
kind of decoherence have been seen in a recent experi- 
ment by Meekhof et. al. 0. However, they expect to sig- 
nificantly reduce decoherence from classical noise sources 
in future experiments |^^ . We thus do not include any 
decoherence effects in our model. 

Treating the laser excitations classically, the two elec- 
tric fields are described by 

^i2(i, t) = ei2 {Ei2 e-^I'^-*— *1 + h.c.} , 

^23{y, t) = {^^23 g-^I'^-^--*! -f- h.C.} , (4) 

where and £23 are polarization vectors, fci2 and k2z 
are wavenumbers, and W12 and W23 are the frequencies of 
the lasers. We assume the laser phases to be absorbed in 
the complex amplitudes E12 and £'23- In dipole approxi- 
mation this leads to the Hamiltonian 

H= ;iwi|l)(l| -t-?iw2|2)(2| -t-?itj3|3)(3| 
) + hi^b{b^ b) 



(a) (b) 
|2) 




FIG. 1. Tlie two-mode Raman transition which couples 
the electronic and the two motional degrees of freedom in the 
X- and y-direction. The effective three-level ion shown in (a) 
is confined within a two-dimensional harmonic trap. As il- 
lustrated in (b) two laser beams propagating in the x- and 
2/-directions generate a stimulated Raman transition between 
the ground states [1) and |3). 



D_12-S.12 i223-£.23i (5) 

where we have denoted the dipole moments of the |1) 
1 2) and |2) <^ |3) transitions by D.12 and 0,23 respec- 
tively. The frequencies uji, lu2 and ^3, are associated with 
the energies of the electronic states |1), |2) and |3), and 
the operators a(a^) and fe(fe^) are the annihilation (cre- 
ation) operators for vibrational quanta in the x- and y- 
directions. These operators are related to the position of 
the ion in the x-y plane through 

X = Axo {a + a^) , 

y = Ayoib + b^), (6) 

where Axq — (?i/2i/am)^/^ and Ayo = (?i/2j/fcm)^/^ are 
the widths of the ground state in the two-dimensional 
harmonic oscillator potential in x- and y-directions, and 
m is the mass of the ion. If the laser beams are sufh- 
ciently far detuned, i.e. 

IA12I, IA23I » |gi2|, I523I, |Ai2 - A23I , (7) 

the two ground states |1) and |3) are coupled via a stim- 
ulated Raman transition and the excited state |2) can be 
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adiabatically eliminated. In the above inequality we have 
defined the laser detunings A12 = (w2 — wi) — a;i2, A23 = 
([^2 — W3) — ^23, and the dipole coupling constants 512 = 
{mi2-ei2\'2)Ei2/n, 523 = {m23-e23\2) E23/n. As dc- 
scribed in Appendix A, the adiabatic elimination proce- 
dure leads to the Hamiltonian 

H= |1>(1| + ?iW3 |3)(3| 

+hva (a^ a) + tivb (p^ b) 

-%i3e~'['=i^*~'=^^*-("i^-'^^^''*l (g) |1)(3| 

-%*3 e'['=i2*-'=23y-(^i2-a>23)t] ^ |3^)^i| (8) 

where we have dropped the term describing the free en- 
ergy of the excited state |2) as in the far detuned limit 
(0) the excited state is no longer connected to the two 
ground states. Furthermore, we have defined the Raman 
coupling constant 



ffl3 = 512523 



1 

All 



1 

A^ 



(9) 



and the energies Tiuji and TicDs of the ground state lev- 
els |1) and 1 3), which are Stark shifted as a result of the 
adiabatic elimination of the excited state, are 



LOl — LOi 



UJ3 = i03 



2|gi2 



^12 



213231 



i23 



(10) 



In order to proceed, we will consider the Raman cou- 
pling Hamiltonian (|8|) in the interaction picture of Ho = 
huji |1) (1| -I- huj^ |3)(3| -I- hva {a^ a) + TiVb (b^ 5), and trans- 
form to the new Hamiltonian 



_ iHot/n 



Hi= e 



{h - Ho) 



^-iHot/h 



(11) 



In doing so and replacing the position operators x and y 
by we obtain the interaction Hamiltonian 



Hi = -?i5i3 exp 



o('7l2 



?723) 



|1>(3| 



(-*r?i2)'"+^ {iV23) 



m! ^! 



exp i {yaim — y\ + Vb[u — n] + t + h.c. , 

(12) 



where we have defined the Raman detuning, 

Ai3 = W12 - ^23 - (^^3 - Wi) , 



(13) 



and the Lamb-Dicke parameters in x- and y-direction, 
7712 = Aa;o/ci2, and 7723 = A7/ofc23- The square of the 
Lamb-Dicke parameter gives the ratio of the single pho- 
ton recoil energy to the energy level spacing in the har- 
monic oscillator potential. 



III. SPECIFIC COUPLING SCHEME 



In this section we construct a particular configuration 
of Raman lasers to decouple the electronic and motional 
dynamics of the trapped ion for suitably-chosen initial 
electronic states. This is done by symmetrically combin- 
ing two Raman transitions as described below. We then 
obtain the Hamiltonian (|l|) in the Lamb-Dicke approx- 
imation and in the limit of suitable trap anisotropy for 
specific sideband detunings of the lasers. 

The electronic and motional dynamics can be decou- 
pled in general for the Hamiltonian i?i = M (g) |1)(3| -I- 
(g) |3)(1|, where M may be any operator that acts 
on the motional degrees of freedom only. This is done 
through the addition of another interaction generated by 
Hj' = M (g |3)(1| -h Aft (g, |1)(3|. Combining both inter- 
actions, we have ~ Hi-\- H[, so that the combined 
Hamiltonian, i7f * = {M + M^) g) (|3)(1| -I- |1)(3|), fac- 
torizes. For the case where Hi is given by (p^), H[ can 
be generated by an extra pair of Raman lasers with suit- 
able detunings, propagation directions and phases. To 
be more specific, we require a symmetric combination of 
two Raman transitions, so that 



K3 


- -Ai3 


V'l2 


to 


V23 


= -^723 , 


913 


= 5l3 7 



(14) 

(15) 
(16) 



where all quantities without primes correspond to the 
first pair of Raman lasers and all primed quantities re- 
fer to the second pair. If, for the first pair of lasers, 
the Raman detuning A13 is given by jl^), then the 
first condition ([l^), requires an appropriate choice of 
the frequencies 07^2 and 1.^23 fo'' the second pair, so that 
A']^3 — uj'12 ~ ^'23 ~ {^^3 ~ '^1) — — A13. This is illustrated 
in Fig.^. The second condition (p^), is satisfied by choos- 
ing the second pair of beams to be counter-propagating 
with respect to the first pair, so that = — fci2, and 



-k23, as seen from the definition of the Lamb- 



Dicke parameters, 7712 = Axofci2, and 7723 = Ai/o^23. 
Here, we have neglected the differences |A;i2| — |^i2li 



and |A;23| 

|t^23 



|fc23|, since \uJ12 



^121 



12; 



and 



■'231 



<C W23,'^23- This restriction can be lifted if 
one chooses the second pair of lasers to be not exactly 
counter-propagating with the first. The third condition 
(|l6| ) , requires a suitable choice of laser phases for the two 
pairs of Raman beams which can be easily read from (|^). 
The symmetric combination of the two Raman transi- 
tions as specified by (p^)-(p6[) then leads to the interac- 
tion Hamiltonian 



3 



= <^ 513 exp 



-0(^712 



^23) 



m! /z! 



I 7y' 



X exp 



i {va[m - fJ.] + VbV' - n] + A13) t 



®{|1)(3| + |3)(1|} 



(17) 



which factorizes. We now assume the ion to be initially in a direct product of its motional and electronic state with the 
electronic state prepared as |+) = (|1) + |3))/V2. This superposition state |+), can be prepared from the ground state 
jl), by applying a resonant 7r/2-pulse (A13 = 0), if the ion is confined within the Lamb-Dicke limit, i.e. 7712,^723 ^ 1, 
1^. The dynamics generated by eg. ([l7|) acting on this state factors and leaves the electronic state unchanged. This 
allows us to reduce the dynamics to that of the motional degrees of freedom only, and we write 



Hl°' = -%i3 exp 



-o('?12 



'?23) 



n ' 



exp i {va[m - /i] + viy[v -n] + A13) t 



h.c. 



(18) 




FIG. 2. Schematic diagram of two symmetric Raman 
transitions which in combination decouple the electronic and 
motional dynamics of the trapped ion for suitably-chosen ini- 
tial electronic states. In (a), the frequencies uju and 1^)23 of 
the two lasers that generate the stimulated Raman transi- 
tion between the ground state levels |1) and |3) are chosen 
such that the Raman detuning A13 = UJ12 — W23 — (io-j, — lui) 
is positive. The coupling lasers are red detuned with 
respect to the |1) <^ |3) transition. In (b), we show 
the symmetric Raman transition to (a). The frequencies 
fi'iid tJ23 of the coupling lasers are adjusted so that 
A'ls = uj{2 —^23 — (w3 — wi) = — A13. The coupling beams are 
blue detuned with respect to the |1) <^ |3) transition. For the 
two transitions to be symmetric we additionally require the 
coupling beams in (b) to be counter-propagating with respect 
to the beams in (a). 



We now discuss the sideband detunings which, in the 
Lamb-Dicke approximation and in the limit of suitable 
trap anisotropy, lead to the desired interaction (|l|). In 
particular, detuning the two pairs of Raman lasers to 
specific vibrational sidebands allows us to choose spe- 
cific values for ka and kh in (|l|). Since we require the 
two Raman transitions to be symmetric, it is sufficient 
to consider the first pair of Raman lasers. Therefore, we 
return to the vibronic Raman coupling Hamiltonian (12). 



From eq.(y_2|) it is clear that by fixing the size of the de- 
tuning A13, i.e. by choosing the frequencies of the two 
coupling lasers, we can tune to a resonance between spe- 
cific vibronic levels. As illustrated in Fig.ljwc introduce a 
virtual level |c) with energy TilJc to help visualize the Ra- 
man transitions between the ion's vibronic levels. If we 
letcji2 = {LOc-^^ij-kaVa-, and [J23 = {^c-^^ij-hvb-, then 
with respect to level |c), the first laser is tuned to the ka- 
th red sideband of the ion's vibration in the a;-direction, 
the second laser is tuned to the fcf,-th red sideband of the 
vibration in the y-direction and the Raman detuning is 



A 



13 



kbVb 



(19) 



This situation is illustrated in Fig.y for the specific ex- 
ample, ka — kb — 1. Now, if only on-resonant terms in 
( p^ ) are retained, we have m = n + ka, and n = + kb, 
and we obtain the Hamiltonian 



= I l)(3|®^%(/z, z.) a^'^'a^V 6^ V 6'^'' + h.c. 



(20) 



where we have defined the coupling constants 

1 



^.913 exp 



:(^?12 +'?23) 



(-i7?i2)^^'+'^'° (^??23)^''+^'' 

fi\{fi + kay. iy\{i^ + kby. 



(21) 



This is a two-mode generalization of the non- 
linear Jaynes-Cumniings model introduced by Vo- 
gel et. al. It is important for the trap fre- 

quencies Va and Vb to be non-commensurate to ar- 
rive at this result. This becomes clear from Fig||. 
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FIG. 3. Schematic diagram of the vibronic energy levels 
that are connected by the two laser beams which generate 
the stimulated Raman transition. The two Raman lasers are 
tuned such that the Raman detuning A13 — kti^t — kai^a, 
with ka = kb = 1. With respect to the virtual level |c), 
the laser propagating in the a:-direction is tuned to the 
first red sideband of the ion's vibration in the 3;-direction 
and the laser that propagates in the j/-direction is tuned 
to the first red sideband of the ion's vibration in the 
y-direction. This causes a resonant transition between the 
vibronic states \na — l)a\'nb)b\'i) \na)a\nb — 1)6|1), where 
the states \na)a\nb)b denote the usual number state basis for 
the two-dimensional harmonic oscillator and the numbers Ua 
and Ub give the number of vibrational excitations in the x- 
and 2/-direction respectively. 



If the trapping potential is isotropic, Va — i^b, and the en- 
ergy levels become degenerate. Consequently, the Raman 
transition Hamiltonian ( p^ ) contains on-resonant terms 
in addition to the ones retained in (|20|). In the exam- 
ple ka — kb = this leads to a coupling Hamiltonian 
Hi cx + {a^ b+ab'< -a'' a-b'' b) + 0{r]^)]<»\l){3\+h.c. , 
where we have assumed the Lamb-Dicke parameters to 
be of the same order of magnitude, 7712 ~ 7723 ~ V- In gen- 
eral, if the frequencies Va and Vf, are commensurate the 
Raman transition Hamiltonian ( ]T^ ) contains resonances 
in addition to the ones considered in eq.(pO|). As we will 
show in the next section, in the Lamb-Dicke limit, the 
coupling constants corresponding to these additional res- 
onances can be greatly reduced by increasing the ratio of 
the trap frequencies, Va/vb- 



The symmetric Raman transition ( |20| ) is generated by 
a second pair of lasers as specified in eqs.(p^)-(p^. In 
particular, we note that ( [l4| ) can be satisfied with the 

choice oj'12 = {'^c — ^i) + kaVa, and W23 = {'-^c — ^3) + kbVb, 
for the frequencies of the second pair of lasers. With re- 
spect to the virtual level |c), these lasers are then detuned 
by the same amount as the first pair but to the blue vi- 



brational sidebands rather than the red. Combining both 
Raman transitions, we obtain the reduced Hamiltonian 

Fft ^ ^t'^-'' fe^V b^b + h.c. , (22) 

for the motional dynamics of the trapped ion as discussed 
above. 

In the last step, we now assume the Lamb-Dicke limit, 
where 7712, 7723 <C 1. In this limit we approximate eq.(|2^) 
by keeping only the lowest order terms in 7712 and 7723. 
From cq. (pi|) , these are the terms /i = = 0, and we 
obtain 



jjtot 



9 a 



b^b 



9 



i^- b^ 



(23) 



where g = (7(0,0), is given in (pi]). The above Hamil- 
tonian (H^) realizes the desired interaction ([|), between 
the two modes a and b of the ion's motion in the x- and 
y-direction. We note that the coupling constant g de- 
pends on the Lamb-Dicke parameters through the factor 
^12 %3- Consequently for fixed laser power, i.e. fixed |gi2| 
and 1.92317 and small Lamb-Dicke parameters, the cou- 
pling strength may be very small. One can increase the 
coupling constant g by increasing the laser power while at 
the same time maintaining inequality (^. This permits 
us to ignore the spontaneous emission from the excited 
state 12) on a timescale 



T ^spont 



\9l2\ 
^12 



\9: 



231 



^23 



(24) 



where 7 is the rate of spontaneous decay from level |2) 
. This is important as the decoupling of the motional 
and electronic dynamics relies on maintaining the coher- 
ence of the electronic degrees of freedom. In section ^ 
we will compare the timescales for spontaneous emission 
and the Raman-generated motional dynamics for the spe- 
cific case of rotation (^), given the parameters of recent 
experiments |7j. 



IV. LIMITATIONS 



In this section we further discuss the approximations 
under which the Hamiltonian (p3|), gives a valid descrip- 
tion of the system dynamics. First, we address the size 
of the corrections that we have neglected in the Lamb- 
Dicke approximation. We then show that the coupling 
constants of the additional resonances in the case of com- 
mensurate trap frequencies can be made as small as these 
corrections for a suitably large ratio of the trap frequen- 
cies, Vajvb- Finally, we discuss the limitations imposed 
on our Hamiltonians from neglecting off-resonant transi- 
tions. 
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A. Lamb-Dicke approximation 



From the previous section it is clear that the Lamb- 
Dicke hmit is an important requirement for us to engineer 
the desired interaction (p3|). The Lamb-Dicke approxi- 
mation led us from (22) to ( p3| ) under the assumption 
?7i2, '723 <C 1. We note that both ( |2^ ) and ( p3| ) couple the 
same vibrational states 



\c) 



b/b 



ka)a\n) 



(25) 



where \m)a\n)b, denotes the usual number state basis for 
the two-dimensional harmonic oscillator. Therefore, we 
do not neglect any additional resonances between other 
states than the ones given in (|2^) by making the Lamb- 
Dicke approximation. 

We define the Lamb-Dicke approximation for suitably 
small 7712, ?723 to be the approximation where all terms 
in (IH) of order 77^ smaller than the leading term are ne- 
glected i.e. 



|g(M,'^)l 
|5(0,0)| 



(26) 



where we have assumed the Lamb-Dicke parameters to 
be of the same order of magnitude, 7712 ~ 7723 ~ "H- 

It is important to note that the orthogonality of the 
Raman laser beams shown in Fig.0 is not essential. In 
fact, the size of the Lamb-Dicke parameters can be 
reduced by changing the geometry of the lasers and 
choosing the two Raman beams to be almost counter- 
propagating. In this situation the wave vectors and 
^23 of the two Raman beams have to be added and the 
numbers ki2 and ^23 in are then the projections of 
k = fci2 + .^23 onto the x- and y-axes respectively. 



B. Trap anisotropy 

As we have mentioned in the previous section, even 
in the case of an anisotropic trap, there are on-resonant 
terms in addition to the ones included in (^0|) when the 
trap frequencies are commensurate. This is illustrated 
in Fig.^ where Va = 5z/b, and again, ka = kh = 1. In 
addition to the |m — l)Q|n)b|3) 4^ \m)a\n~ l)h|l) tran- 
sition shown in Fig.|[ the \m)a\n ~ 4:)b\<i)<^ I'm) a\n) b\l) 
transition is resonantly coupled as in FigJj. In the follow- 
ing we show that in the Lamb-Dicke limit, the coupling 
constants g, corresponding to these additional resonances 
satisfy 



\9\ 



|5(0,0)| 



< 0(7y2) , 



(27) 



if the ratio of the trap frequencies is chosen 
large enough. These additional terms can thus 
be neglected in the Lamb-Dicke approximation. 




FIG. 4. Vibronic energy level diagram for the case of 
the two trap frequencies Va and being multiples of each 
other, Va = 5ut. As in Fig| the two Raman lasers are 
tuned such that the Raman detuning A13 = ktVt — kaVa, 
with fca = fci, = 1. In addition to the desired resonant tran- 
sition \na ~ l)a|ni,)6|3) <^ \na)a\ni, - 1)6|1>, shown in grey, 
the \na)a\nb — 4)i,|3) <^ |wa)a|n6)6|l) transition is resonantly 
coupled as shown in black. In the Lamb-Dicke limit, the cou- 
phng constant corresponding to this additional resonance can 
be reduced to the size of the corrections to the Lamb-Dicke 
approximation for the desired resonance by increasing the ra- 
tio of the trap frequencies, Va/fb- 

We start by deriving the resonances that occur if the 
two trap frequencies Va and Vf, are multiples of each other. 
Without loss of generality we choose 



= lVb; 



(28) 



where / is a positive integer number. In deriving the in- 
teraction (^^ the laser frequencies were chosen to give 
the detunings A13 = k^Vb - kaVa, and A'^g = kaVa - k^Vb, 
for the two pairs of coupling beams respectively. We will 
explicitly consider only the first of these two cases, i.e. 
A13 = feftj/fc — kai'a, since the second follows analogously 
by interchanging the operators a <^ , and 6 <^4> 6^ , and 
leads to the same limits for the trap ratio I = Va/vb- Now, 
with (p8|), the resonance condition in ( [l^ ) becomes 



(n 



0, 



(29) 



where all numbers are positive integers. In order to sim- 
plify the discussion we categorize the resonances by in- 
troducing an integer number N, and rewrite (E9h, so that 



II — m = —ka + N , 
h' — n = —kii + I N . 



(30) 



Following this categorization, we divide the resonances 
determined by (p9|) into the three cases, [1] N — 0, [2] 
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N positive, and [3] N negative. We subdivide case [2] 
further into 

[2a] (i) 0< N <ka and < N I < kt , 
{a) 0< N <ka and N I > h , 

[2b] {i) N >ka and < N I < k^ , 
{a) N > ka and N I > k^ . 

Below we will examine the cases [1] and [2a]{n) in detail 
as the later case contains resonances with the largest con- 
tribution besides the required resonance at iV = 0. We 
have examined the other cases and will not repeat their 
analysis except to note that they all give rise to leading 
order corrections of order higher than those found in case 
[2^{ii) in r]. Thus, to obtain the desired Hamiltonian 
(P3|), the resonances in case [^a](ii) will impose the most 
stringent condition on the size of the trap ratio / = Vajv^,. 
Throughout this discussion we will consider only the low- 
est order terms in the Lamb-Dicke parameters, since we 
have already addressed the size of the corrections to the 
Lamb-Dicke approximation in the above. 

In case [J], we have iV = 0, so that from the reso- 
nance condition in ( |30| ) we obtain n — m — ~ka, and 
V — n = —kf,. This is the case which leads us to the de- 
sired interaction (|2^), and which we have discussed in the 
previous section. 

We now consider case [2a]{ii). Here, we have < < 
fca, and N I > kb, so that from (pO|) we obtain 



'm = fi + ka — N>fi, 
i> = n + lN — ki,>n. 



(31) 



Inserting these identities into (^^ and keeping only the 
lowest order terms in the Lamb-Dicke approximation, i.e. 
/Lt = n = 0, we obtain 



where we have defined the coupling constant 



(32) 



-513 exp 



\lN-ku 



{ka-N)\ {lN-kb)\ ■ 
(33) 



We require the coupling constants of the above reso- 
nances ( |3l| ) to be smaller than or equal to the coupling 
constants of the terms that we have neglected in the 
Lamb-Dicke limit (1271). Therefore, we have the condition 



\9\ 



knl 



k,\ 



|5(0,0)| {ka^N)\{kt-lN)\ (7712) 



/„ \lN-2ki, 
{V23) < ^2 



N 



(34) 



where again, we assume both Lamb-Dicke parameters to 
be of the same order of magnitude, 7712 « 7723 ~ rj- In 
order to derive a limit for the trap ratio I from the above 



expression, we consider the factor containing the Lamb- 
Dicke parameters and the one containing the factorials 
separately. If r/12 ~ r/23 w 77, we have 



('723) 



lN~2ku 



{ni2) 



N 



N{l-l)-2k^ < jj2 ^ 



(35) 



which is satisfied if, N{1 — 1) — 2fcb > 2. Since this condi- 
tion has to hold for all N in the range < iV < fc^, this 
leads to the requirement 



/ > 2A:h -I- 3 , 



(36) 



for the trap ratio I. Next, we consider the term including 
the factorials. We require this term to be smaller or equal 
to unity as under the above condition (^) the factor con- 
taining the Lamb-Dicke parameters already satisfies (|3^). 
We have 



khl 



< 



kj.ktl 



< 



{ka + fcb)! 



{ka-N)i {kb-iNy. - {iN-kby. - {IN -kbY. ' 

(37) 

where in the first inequality we made use of the fact that 
for the resonances we are discussing here, < N < ka, 
and the second inequality holds since, {ka + kbY- > kal kbl, 
for all positive integers ka and kb- From the above in- 
equality (^7|) the factor containing the Lamb-Dicke pa- 
rameters is smaller or equal to unity if, {I N ~ kbY- > 
{ka + kbY, ^^'^ since this has to be satisfied for all N in 
the range < TV < ka, we require 



l>2kb + ka. 



(38) 



Depending on the interaction that we want to generate, 
i.e. depending on the number ka, the inequality ( p6[ ) or 
( p8| ) will impose the stronger limit on the trap anisotropy. 
For the two examples given in (|^) and (|), we have 
ka = kb = 1, and ka = 3, kb = 1, respectively. Therefore, 
in order to generate the linear coupling Hamiltonian (|^) 
we require the trap ratio I = Va/^b ^ 5, (|3^). For the 
cubic interaction (^) a trap ratio of I — Vajvb > 5, is 
needed from (^. 

For the remaining cases [2a]{i), [2b]{i), [2b]{n) and [3] 
a similar analysis shows that the requirements (36) and 



( pq ) are sufficient to limit the strength of these resonances 

Although in the above discussion we have explicitly 
assumed the two trap frequencies Va and i>b to be mul- 
tiples of each other, the limits (^6|) and ( |38| ) also hold 
for commensurate trap frequencies. In this case the trap 
ratio is a rational number, i.e. I — p/q, where p and q are 
positive integers. Since in the resonance condition ( |30| ) 
all numbers need to be integers, the number N which 
categorizes the resonances can only take on multiple val- 
ues of q, so that IN = pN/q, is an integer. As we have 
discussed all integer values of N, any trap ratio I = p/q, 
which satisfies the inequalities ( |36| ) and (|3^) suffices for 
the unwanted resonances to satisfy (p^). Hence, for given 
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values of ka and kb, the coupling constants of all addi- 
tional resonances due to energy level degeneracies in the 
case of commensurate trap frequencies are at least a fac- 
tor of 77^ smaller than the coupling constant of the desired 
resonance (p3|), if the trap ratio is chosen large enough 
according to the limits in (|3^) and (^sj). 



C. Off-resonant terms 



r|5(o,o)p 



MINI 



{M - kaV. (N - kb + ly. \m3 



kh 



— (40) 



where we have assumed the Lamb-Dicke limit to calculate 
the coupling V between the states (|39|). We will further 
investigate the significance of the limitations discussed 
here in the section below where we concentrate on the 
linear coupling Hamiltonian given in (^). 



As pointed out by Gardiner et. al. ||25|, dropping all off- 
resonant terms in going from ( p^ ) to (|23|) iinposes a limit 
on the time T for which the Hamiltonian (|2^) is a valid 
approximation. This limit can be calculated in second 
order perturbation theory to be T F^/A ^ 1, where V 
is the effective coupling to the nearest off-resonant tran- 
sition in (|8|) and A is the corresponding detuning. If 
\M)a\N)b is a characteristic state which represents the 
highest energy state that we allow to be acted upon, the 
transitions 

\M-ka + l)a\N)b^\M)a\N~kb)b, 

\M ~ ka)a\N)b ^ \M)a\N ~ kb + l)b (39) 

are the strongest coupled off-resonant terms. For these 
two transitions the limit becomes 



7^15(0,0)1^ 



Ml m 



{M - ka + ly. (N - kbY. \m2 



V. ENGINEERING ROTATION 



In the following we use the above formalism to tar- 
get the linear coupling Hamiltonian Hi^\ given in (|^) 
and show how this generates a rotation of the two- 
dimensional quantum motional state of the ion. We then 
examine the validity of the approximations discussed in 
the previous section through a numerical analysis of this 
specific example. 

The linear coupling H^^^ is obtained from the sym- 
metrically combined two-mode Raman Hamiltonian jl^ ) 
through the particular choice A13 = Ub — Va, for the Ra- 
man detuning and adjusting the relative phase of the 
lasers such that the Raman coupling constant 1713 = 
1 1 (713 1, is purely imaginary. This leads to the Hamilto- 



Hl°' ^ -ih\gi3\ exp 



^-^ to! fj,\ n\ v'. 



X exp 



« (i/q [m - ^ - 1] + i^b [i^ + 1 - n] + A13) i 



h.c. 



(41) 



which, in the limits discussed in the previous section, re- 
sults in the linear coupling h[^\ The coupling constant 
g in (^) is then given by 17 = — [(/(O, 0)|. The Hamilto- 
nian iJj effects a rotation of the two-dimensional quan- 
tum motional state of the ion about the centre of the 
trap. This can be seen by examining the action of the 
Hamiltonian h\^^ on the operators a and h. Using the 
Baker-Campbell-Hausdorff theorem, we have 

ag = C/(i)at/(^'^ =acos6l-6sin6l, 

&e = C/'^^foC/^^^^ = asin6l + &cos6', (42) 

where the angle 6 = gt, and C/*-^^ is the unitary transfor- 
mation generated by the Hamiltonian h\^\ i.e. 



(43) 



From mj it is clear that the transformation in (E2) cor- 



responds to a rotation of the rescaled coordinate system 
X = x/ Axq, and y = y/Aj/o, through an angle 9 = gt, so 
that in the rotated coordinate system we have 



xe = X cos u — y sm c 
yg = xsvaO + ycost 



(44) 



Now an arbitrary pure or mixed motional state of the 
ion is characterized by a density operator p which can be 
written as 



P'^'^n \m)a\n)b{lJ-\a(l^\b 



/ J rrn.n 



V ml nl 



|0)a|0)6(0|a(0|6 



(45) 



The time evolution of this state under the action of the 
Hamiltonian H^^^ is then given by 
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E 



rm.n 



4.m "4.71 „ - 

an bl a^6^ 



-\/m! 7i! 



|0)a|0)b(OU(0|t 



6 J 



(46) 



where we have used (^2|) and |?7T-)a|?i)b is the number 
state basis for the two-dimensional harmonic osciUator 
but now in the rotated coordinates xg and yg as given 
in (^^. Therefore the motional state of the ion given 
by p{t) is identical to p, but rotated through an angle 
9 — gt. In particular, this is accomplished without prior 
knowledge of the motional state p. 

Having convinced ourselves that the linear coupling 
Hamiltonian H^"^ does rotate an arbitrary motional state 
of the ion we now examine the validity of the approxi- 
mations discussed in the previous section for this specific 
example of the general coupling Hamiltonian (^3|). We 
consider a state rotation through the angle 6 = 7r/2, so 
that Tj-ot = 7r/2(/ is the required time to rotate the state. 
For this case , the limitations due to off- resonant terms 
(section IV C| ) as given in eq.(40) take the form 



2 I5i3| 



(47) 



in the limit of small 7712 ~ 7723 ~ Jy, and where we have 
assumed to be the smaller of the two trap frequencies. 
Here, iV^ax = max(Af,M) in (|^). From eq.(|^ it is 
clear that the ratio of the lower trap frequency over the 
Raman coupling constant, 



7 = Vbl\9i^\ , 



(48) 



determines the significance of off-resonant terms in the 
system dynamics. 

From our discussion of the si gnificance of additional 
on-resonant terms (section IV B ) we require a trap ratio 
I — Va/vb > 5, for the linear coupler where ka = kb = 1, 
(^) . The estimates used to determine this minimal trap 
ratio essentially compare the coupling strengths of the 
different terms appearing in the Hamiltonian ( [4l| ) with no 
reference to the actual state on which it acts. Although 
this method of estimation is used in the literature, it can 
only serve as a rough guide. A more rigorous measure 
of how the unitary time evolution ^7^^^ , generated by the 
symmetrically combined two-mode Raman Hamiltonian 
(p]), deviates from the desired unitary evolution 
generated by the linear coupling Hamiltonian (|2[), can 
be quite complicated. A fully rigorous state independent 
measure of the difference between two unitary operators 
can be constructed but we will not consider this 

here. In order to examine the validity of the approxi- 
mations discussed in the previous section we adopt the 
overlap 



j^Kvi/totl^-)!, 



(49) 



as a measure of the deviation between the two unitary 
evolutions U^ot ^'^'^ U^^^ for an initially pure quantum 
state l^'o)- Here the state 



(50) 



gives the unitary evolution of the initial state |^'o) under 
the action of the symmetrically combined Raman Hamil- 
tonian (p4|), and the state 



I*; 



0; 



(51) 



gives the desired evolution of the initial state under 
the action of the linear coupling Hamiltonian (^). This 
cannot be calculated analytically. To go beyond the ana- 
lytics we numerically compute the unitary evolution ([5^), 
including the higher on-resonant and off-resonant terms, 
on the initial pure state, = \a)a ® |a)b, where \a)a 
and I a) a are coherent states in the vibrational modes a 
and b respectively. In this case the desired state (|5l]), 
rotated through 6 — 7r/2, is given by \^) = \—a)a ® \a)b- 
The results of our numerical analysis are shown in FigJq. 




FIG. 5. Results from our numerical analysis of the devia- 
tion of the unitary evolution generated by the symmetrically 
combined Raman Hamiltonian tuned for rotation from the 
desired state rotation. We plot the overlap 5 = |(^tot|^')|, be- 
tween the state l^tot), resulting from the Raman Hamiltonian 
tuned for rotation and the desired state l^*) = | — a)a ® \ct)b, 
resulting from a rotation of the initial state j^o) = \ct)a®\ct)b, 
through the angle 6 = 7r/2. We have chosen a = 1. The graphs 
(a) - (d) show the dependence of the time evolution of 5 on 
the parameter 7 which takes on the values 7 = 22/2" where 



There we plot the overlap (5 as a function of the scaled 
time gt, for different values of the parameter 7, given 
by ( ^8| ) and a coherent state amplitude a = 1. Be- 
fore we discuss our results we note the following on our 
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choice of parameters. In Fig.||(a) - (d), 7 takes on the 
values 7 = 22/2", where n = 3,..,0 in unit steps, and 
the Raman couphng constant 1713 is kept constant. For 
simphcity, we assume the geometry of the laser excita- 
tion to be arranged so that the Lamb-Dicke parameters 
7712 = Axofci2, and 7723 = ^Uokis, are equal. Here it 
is important to note that the values of the Lamb-Dicke 
parameters 7712 and 7723 depend on the size of the trap 
frequencies Va and Vb through Axq = {h/2vamy/^ and 
Aj/o = (?i/2z/f,m)^/^. Therefore the size of the Lamb- 
Dicke parameters depends on 7 and varies from Fig|^(a) 
to (d). To incorporate this dependency in our numer- 
ical analysis we put 77;^2 = IV23 — 0.88, which gives 
'712 = V23 = 0.2, when 7 = 22. These are values for 
the Lamb-Dicke parameters and the ratio 7 which have 
been demonstrated in cold ion experiments 0-^1 . Fol- 
lowing our discussion of the trap anisotropy we choose 
the trap ratio / = i^a/i^b = 5. Our numerical analysis 
was performed in a finite (truncated) number state ba- 
sis (|0)a|0)f,..|8)a|8){,) with a cutoff chosen such that an 
increase of this cutoff does not significantly alter the re- 
sult of our integration. Fig.^(a) shows the time evolution 
of the overlap 5 for the lowest value of 7. Here the off- 
resonant terms in eq.(pl|) cause strong modulations in 6. 
For higher values of 7 these modulations become much 
less pronounced as the off-resonant terms contribute less 
on these timescales (Fig.||(b) - (d)). The time evolution 
of 5 for the highest value of 7 is shown with the solid line 
in Fig.||(d). The plot reaches a maximum oi 5 ~ 0.99 at 
gt « 1.02 X 7r/2. It shows almost no deviation from the 
dashed line in Fig.^j (d) which is a numerical integration 
of the system dynamics where we only include the desired 
resonances in the Hamiltonian as in (^2|) . The numerical 
analysis shows that for low system excitation the Hamil- 
tonian can be engineered with high accuracy within 
present ion traps. To investigate the system dynamics 
for higher energies becomes computationally very expen- 
sive. To achieve the same accuracy as obtained for a = 1 
for higher values of a, the number state basis must be 
greatly enlarged. 

In the above, we neglected decoherence. From our final 
comments in section [II, using the experimental parame- 
ters M for 9Be+ with 7/27r = 19MHz, 5i3/27r = 500kHz, 



Ai2/27r = A23/27r = 12GHz and 7712 = 7723 = 0.2, we find 
Tspont ~ 200/is. This is to be compared with the time to 
rotate the motional state through the angle 6 = tt/2, 
Trot ~ 12fis. This confirms our initial assumption that 
decoherence through spontaneous emission can be ne- 
glected for this process. However, this may not be the 
case when engineering higher order interactions. One can 



shorten the interaction time by increasing the laser power 
while maintaining inequality (|^) by increasing the detun- 
ings A12 and A23. The fundamental limit is then given 
by the detunings that one can realize and the accessible 
laser power. 



VI. CONCLUSION 

In this work we showed how one can engineer a class of 
Hamiltonians for the motional dynamics of an ultra-cold 
ion in a harmonic trap. The process uses a stimulated 
Raman transition in a A-configuration with the two lasers 
propagating along the x- and y-directions. To decouple 
the internal electronic dynamics from the external mo- 
tional dynamics we constructed a Hamiltonian in which 
these evolutions factored. This was done through the 
addition of a second pair of lasers which generated the 
symmetric counter-part to the Hamiltonian generated by 
the first pair of lasers. By preparing the electronic states 
in a particular superposition, the internal and external 
dynamics completely separated and we could treat the 
motional dynamics alone. In the Lamb-Dicke limit and 
with suitable sideband detunings, we could "target" a 
particular term to be of leading order in the Hamilto- 
nian. However, we found that besides the term we wanted 
to dominate, other, higher on-resonant terms appeared. 
We could manipulate the strengths of the couplings to 
these unwanted terms by altering the trap frequency ra- 
tio and found that we could neglect these unwanted terms 
in the Lamb-Dicke approximation for large enough trap 
anisotropics. Finally, we did a numerical evaluation of 
the full Hamiltonian as a check on the analytical esti- 
mates. Although we have primarily concentrated on the 
linear rotation Hamiltonian (Q), higher order dynamics 

can be generated i.e. H^^^ , given by eq. (||) . The nonlinear 

" (3) 

Hamiltonian , has been much studied in the quantum 
optical literature as a model of nonlinearly-coupled field 
modes |2^] . We know from this work that such Hamilto- 
nians generate a rich nonlinear dynamical structure re- 
fiecting the strong mode entanglement characteristic of 
those couplings. Their optical realization is difficult, but 
may well be more straightforward in trapped ion dynam- 
ics, as resonances can then be used to isolate chosen non- 
linearities. 

Finally, we note the recent publication of two papers 
P?! which examine types of non-linear interaction Hamil- 
tonians in the motion of trapped ions which are closely 
related to the work presented here. 
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APPENDIX A 

In this appendix we derive the effective Hamiltonian given in (^) which follows from the adiabatic elimination of the 
excited level |2) and describes the Raman coupling between the two ground state levels |1) and |3). After performing 
the rotating wave approximation the Hamiltonian in (|^) becomes 

H= + ?iw2|2)(2| + ?iw3|3)(3| + nva{a) a) + hvb{b'^b) 

- |1)(2| ® %i2 e-*('=i=*-'^i^*) - |2)(1| ® hgl^ gi{ki2£~u.i2t) 

- |3) (2| (g) %23 e-' ('=^36-c^23t) _ |2) (3| ® hg^^ (^^^S-^^st) ^ (52) 

where we have defined the dipole coupling constants .gi2 = (l|Z2i2-£i2|2) Ei2/h and 1723 — (3|i223-623|2) E23/h. In order 
to compare the tiniescales of the transitions induced by the two laser beams we consider the Heisenberg equations of 
motion for the transition operators ai2 = |1)(2| and ai^ = |1)(3|. 

i ^ ^12 = (^2 - c^i) ^12 - 9*12 C^-*— (^n - ^22) - 523 C^^^^""^^*) #13 , 

i ^ ^13 - {UJS - C^i) ^13 - 523 ('=^3t-23t) + g*^ {k^.^-^^.t) (53) 

Here, all operators (denoted by overbars) are taken in the Heisenberg picture, i.e. CT12 = U{t)(Ti2W {t), where U{t) = 
Texp [i (J* H{t')dt' 1%) ] is the time ordered evolution operator. Using the transformation 

^12 - e-*"-*?i2, 

^23 = e^"^^*CT23, 

^13 = ^12^23 = e-' (-1^— 23)t ^ (54) 
to remove the explicit time dependencies from ( p3|) we have 

?i2 = A12 5^12 - 5i*2 e^'"^ (^11 - ^22) - 52*3 e*'^-^ ai3 , 

?i3 = (A12 - A23) ai3 - 523 e-''==^^5^i2 + 51*2 e^''^=^5^23 ■ (55) 

Under the assumption of large detunings, as given in (|^), we obtain the adiabatic solution for ai2 by setting d'(Ti2/dt = 
0, 1^ so that after restoring the rapidly oscillating time dependence, we obtain 

^ i^*'^ (fe-*— *) - #22) + .g23 e' ('^-^^ *) #13} . (56) 
For the |2) |3) transition we find in an analogous manner 

^ i {^^23 (^33 - ^22) + 5^2 f"-^-"-*) 5^31 } . (57) 

Upon inserting these adiabatic solutions for a" 12 and CT32 into (^2|), we have 

H = hui\l){l\ + hu:i |3)(3| + hva (at a) + hvb (b^b) 

-hgi3 e"'(''^i^*-'=^^^-('^i^-'^^3)*' (g, |1)(3| - %i*3 e*(fei2*-fc236-(c.i2-t^23)t) ^ |3)(i| ^ (53) 

where we have dropped the term describing the free energy of the excited state |2) since in this adiabatic approximation 
it is no longer connected to the two ground states. Furthermore we have defined the Raman coupling constant as 
given in (||), and the energies huji and huj^ (|o|) of the ground state levels |1) and |3), which are Stark shifted as a 
result of the adiabatic elimination of the excited state. 
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